Hydrodynamics within the Electric Double Layer on slipping surfaces 
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We show, using extensive Molecular Dynamics simulations, that the dynamics of the electric dou- 
ble layer (EDL) is very much dependent on the wettability of the charged surface on which the 
EDL develops. For a wetting surface, the dynamics, characterized by the so-called Zeta potential, 
is mainly controlled by the electric properties of the surface, and our work provides a clear inter- 
pretation for the traditionally introduced immobile Stern layer. In contrast, the immobile layer 
disappears for non-wetting surfaces and the Zeta potential deduced from electrokinetic effects is 
considerably amplified by the existence of a slippage at the solid substrate. 

PACS numbers: 68.15-|-e,47.45.Gx,82.45.-h 



The electric double layer (EDL) is a central concept in 
the understanding of the static and dynamical properties 
of charged colloidal systems. This notion was introduced 
in the early works of Gouy, Debye and Hiickel to de- 
scribe the distribution of micro-ions close to a charged 
colloidal surface. The EDL width determines the electric 
interaction range between macromolecules and therefore 
controls the static phase behaviour of these systems. On 
the other hand, at the dynamical level the EDL is at 
the origin of numerous electrokinetic effects Q : elec- 
trophoresis, electro-osmosis, streaming current or poten- 
tial, etc. Because these various phenomena are governed 
by the surface of the sample via the EDL, they provide 
smart and particularly efficient ways to drive or manip- 
ulate flows in micro-fluidic devices 0, Q , where surface 
effects become predominant. 

The extension of the EDL is typically on the order of 
a few nanometers and electrokinetic phenomena there- 
fore probe the nanorheology of the solvent -|-ions system 
at the charged surface. This raises therefore some doubts 
about the validity of continuum approaches to describe 
the dynamics at such scales. These doubts are partic- 
ularly relevant concerning the traditional description of 
the EDL dynamics, which relics both on the mean-held 
Poisson-Boltzmann theory of the micro-ion clouds, but 
also on continuum hydrodynamics for the flow fields Q. 
These two aspects are embodied in the so-called Zeta po- 
tential, denoted which is traditionally defined as the 
electric potential ^^(^s) computed at the surface of shear 
Zs, where the fluid velocity vanishes 0,11,110. This def- 
inition of however relies on the somewhat uncontrolled 
assumption of no-slip boundary condition of the solvent 
at the solid surface 

This assumption has been critically revisited in the 
last years. Indeed, a lot of progress has been made re- 
cently in the understanding of the rheology of fluids at 
small scales, thanks in particular to computer simula- 
tions, such as Molecular Dynamics (see eg and ref- 
erences therein), but mainly to the development of new 
experimental techniques, such as Surface Force Appara- 



tus (SFA) or Atomic Force Microscope (AFM) H E3- 
The conclusions emerging from these studies are that, 
while continuum hydrodynamics are found surprisingly 
to remain valid up to very small length scales, the no- 
slip boundary condition (BC) for the fluid velocity at 
the solid surface may be violated in many situations (see 
and refs. therein). Moreover, this violation 
of the usual no-slip BC is found to be controlled by the 
wetting properties of the fluid on the solid surface : while 
the no-slip BC is fulhlled on hydrophilic surfaces, a finite 
velocity slip is measured on hydrophobic surfaces j^, (l^l . 

In this Letter we show using Molecular Dynamics 
(MD) simulations that a finite slip effect for the solvent 
at a charged surface considerably enhances the measured 
electrokinetic effects. This results in an enhanced C po- 
tential, the origin of which lies in the dynamics of the sol- 
vent at the surface: the flow-induced current, driven by 
the convection of micro-ions in the EDL, is much larger 
in the presence of velocity slip at the surface than for 
the corresponding situation with a no-slip BC. More- 
over, while an immobile Stern layer develops close to 
the charged surface on hydrophilic samples, and does not 
contribute to electrokinetic effects, this layer is absent on 
hydrophobic surfaces and does contribute in a significant 
way to charge transport. 

We first precise our microscopic model and some de- 
tails of the simulation procedure. The fluid system (sol- 
vent and micro-ions) is confined between two parallel 
solid substrates, composed of individual atoms fixed on 
a fee lattice. The solvent and solid substrate particles 
interact via Lennard- Jones (LJ) potentials. 
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with identical interaction energies e and molecular di- 
ameters a. The tuning parameters allow to adjust 
the wetting properties of the fluid on the substrate |3| : 
for a given fluid-fluid cohesion cpp, the substrate ex- 
hibits a "hydrophilic" behavior for large fluid-solid co- 
hesivity, cfs, and a "hydrophobic" behavior for small 
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cps- Here, the wetting (respectively non- wetting) situa- 
tion is typically achieved by taking cps = 1 (resp. 0.5) 
for a fixed cpp = 1.2. This leads to a contact angle 
6* of a liquid droplet on the substrate, measured in the 
simulations equal to 80° (resp. 140°) for a temperature 
ksT/t = 1 (we refer to Q for an exhaustive discussion 
on this point). On the other hand, micro-ions interact 
both through LJ potentials, as described in Eq. and 
Coulomb potential in a medium with dielectric permit- 
tivity Cd [vapii") = kBTqaqpiB/r where qa, qp are the va- 
lences of the interacting charges, and £b = / {AirtdkBT) 
is the Bjerrum length, e denoting the elementary charge] . 
In water at room temperature Ib = 0.7nm. We shall 
choose in our simulations Ib = o". Ions also interact di- 
rectly with the surface atoms via the same LJ potential as 
the solvent {i.e. same cfs). Micro-ions and solvent par- 
ticles have the same size a (in Eq. (^). Wall atoms are 
organized into five layers of a fee solid (100 direction) in 
both walls. For each wall, the first layer only, in contact 
with the fluid, is charged. The corresponding N.u, atoms 
bare a discrete charge, with valency qwaii = —Z/N^u so 
that each wall bears a negative net charge — Ze. The 
solvent contains Z monovalent counter-ions, to which 
Ns — + salt ions are added, all with valence 
one. Global electroneutrality is enforced by imposing 
N-^- = iV_ . The systems simulated are generally made up 
of 10"* atoms. A typical solvent density is Pfcr^ ^ 1, while 
the concentration of micro- ions ps = N±/V is varied be- 
tween pscr^ = 5.10"'^ and psa^ = 0.16 (with V the total 
volume of the sample). This corresponds typically to an 
ionic strength varying between 10~^M and IM. The cor- 
responding Debye screening length (see below) is accord- 
ingly of the order of a few Bjerrum lengths. The charge 
per unit surface on the wall is S = — 0.2e.(T^^, associated 
to a Gouy-Chapman length A ~ {2ttIb^)^^ — 0.8<j. Us- 
ing a typical value a = 5 A, this corresponds roughly to 
— 0.13C.m~^. Periodic boundary conditions are applied 
in the x and y directions, with = Ly = 16cr while 
the distance between the walls is Lz ~ 20.8ct. Ewald 
sums are used to compute Coulombic interactions (as- 
suming a periodicity in the z direction with a box size 
of 112(7, much larger than the wall to wall distance) [lll |. 
Lennard-Jones units are used in the following [distance 
a, time r — (mcr^/e)^/^]. Temperature was kept con- 
stant by applying a Hoover thermostat to the y degrees 
of freedom, i.e. in the direction perpendicular to the flow 
and confinement. 

Our model therefore includes the discrete nature of 
the solvent and charges and a tuning wettability of the 
surface, whereas these effects are usually neglected in 
the traditional description of electrokinetic phenomena. 
Note that we chose to describe the charge interaction 
at the level of an effective dielectric media (with dielec- 
tric permittivity e^). This simplifying assumption -that 
could be relaxed using a more realistic model for the sol- 
vent 5]- allows to investigate specifically the generic in- 
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FIG. 1: Micro-ionic density profiles, averaged over the xy 
directions {pscr^ = 0.07, wetting case). Symbols : MD simu- 
lations results for the counter-ions (o) and co-ions (□); Solid 
and dashed lines : corresponding predictions of the modified 
PB description (see text). Inset: Electrostatic potential. The 
wall is located at Zu, = —10. So". Symbols (o): MD simulation 
results calculated from Poisson's equation and micro-ions pro- 
files; Dashed line : bare PB prediction (see text). 



terplay between slip effects and electric behavior, which is 
the main focus of this work. We have moreover explored 
situations beyond the ones presented here (with or with- 
out salt, using different cps for the ions and the solvent, 
and various Bjerrum lengths Ib), leaving the present con- 
clusions unaffected. 

We now turn to the simulation results. We first fo- 
cus on the equilibrium properties of the EDL. In Fig. 
n we show typical density profiles. The latter are 
found to exhibit important structuration effects close to 
the charged surface and thus depart strongly from the 
Poisson-Boltzmann (PB) prediction 0. However, the os- 
cillations in the micro-ions profiles originate in the struc- 
turation in the solvent itself, and such an effect can be 
captured by a modified PB description. Indeed, due to 
the presence of the solvent, micro- ions not only orga- 
nize due to electric interactions (which corresponds to 
the usual PB description) but also due to the effective 
external field associated with the structuration in the 
solvent defined as Vext{z) = -fesTlog [p/(z)/p/], 
with Pf{z) the solvent density profile and pf its bulk 
value. The micro-ions density profiles p±{z) correspond- 
ingly obey a modified Boltzmann equilibrium : 

p±(z) OC e'5(TeV(.)-V..*(.)) ^ p^{^^)ef0-V(z) (2) 

with (3 — I/UbT and V{z) the electrostatic potential. 
Actually, such a relationship emerges naturally from a 
simple Density Functional Theory, in which the discrete 
nature of both solvent and charged atoms is taken into 
account exactly, while the standard mean-field PB free 
energy is assumed for the electrostatic part. Using Pois- 
son's equation, the electrostatic potential is found to 
obey a modified PB equation, /3eAV = k'^^{z) smh{(3eV) 
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where = SttIbP^^^^ is the Debye screening factor de- 
fined in terms of the bulk micro-ion concentration, and 
7(2:) = pf{z)/ pf is the normahzed sofoen< density profile. 
In order to test the predictions of this approach, we have 
measured the fluid density profiles, Pf{z), and solved 
Poisson's equation with the microionic densities given by 
Eq. using Neumann boundary conditions, and as- 
suming a smeared (uniform) surface charge on the wall. 
As shown in Fig. ^ this approach leads to results in very 
good agreement with the simulations profiles. Moreover, 
a further approximation can be proposed: The solution 
of the modified PB equation for electrostatic potential is 
actually very well approximated by the "bare" PB solu- 
tion Vpb{z) (corresponding to 7(2;) = 1), whose analytic 
expression can be found in the literature Q. This leads 
to p±{z) oc pf[z)exp[^l3eVpB{z)]. The validity of this 
approximation -surprising in view of the strong layering 
effect at work- is emphasized in Fig. ^ (inset), where the 
corresponding bare PB potential |l| is plotted against the 
"exact" electric potential. The latter is obtained from the 
simulations using Poisson's equation by integrating twice 
the charge density profile pc = &{p+ ~ P-)- 

We now investigate the dynamical properties of the 
EDL. We first consider a streaming current experiment: 
an external volume force, /o, is applied to the fluid in the 
X direction, enforcing a Poiseuille flow in the cell, and 
the electric current, /e, associated with the convective 
motion of the micro-ions is measured. The standard EDL 
description of this clcctrokinetic effect predicts a linear 
relationship between the current and the force, in the 
form : 

le = -—Ah (3) 

■n 

where 77 is the shear viscosity of the fluid and A the fluid 
slab cross area. In the simulations, a force per parti- 
cle is applied to all fluid particles and the corresponding 
electric current is measured. We emphasize that linear 
response (in the applied force) of the system was carefully 
checked. In the following we shall use this expression as 
the definition of the - apparent - C, potential^ in line with 
experimental procedures. 

We first discuss the measured velocity profiles. The 
situation corresponding to a wetting substrate is shown 
in the main plot of Fig. |21(here for /g = 0.02 in LJ units). 
The velocity profile is found to exhibit a parabolic shape 
as predicted by continuum hydrodynamics, even at the 
scale of the EDL. Moreover the viscosity, deduced from 
the curvature of the parabolic shape, is measured to keep 
its bulk value. Nevertheless, the no-slip BC is found to 
apply inside the liquid, at a distance of about one layer 
of solvent particle, in agreement with previous theoretical 
predictions j^]. This position of the no-slip BC here de- 
fines the "plane of shear" position, z^, usually introduced 
in the electrokinetic literature 1] . As shown in Fig. |21 the 
layer of micro- ions located within Zs , does not contribute 



to the convective transport, thereby reducing the global 
streaming current. This immobile layer coincides with 
the so-called Stern layer of immobile micro-ions close to 
the charged surface Note moreover that we found 
Zs to vary slightly with electrostatic parameters (surface 
charge S, Bjerrum length Ib), as expected. 

On the other hand, the non-wetting case exhibits a 
very different behaviour, as shown in the inset of Fig. |21 
First, concerning the velocity profile, a large amount of 
slip is found at the wall surface, in agreement with ob- 
servations on uncharged non- wetting surfaces Q. More 
quantitatively, slippage is characterized by a slip length, 
h, defined as the distance at which the linear extrapo- 
lation of the velocity profile vanishes. This amounts to 
replacing the no-slip BC by a partial slip BC, defined as 
^ff = w at the wall position [sj. As shown in the inset 
of Fig. 12 the velocity profile is found to be well fitted 
by the continuum hydrodynamics (parabolic) prediction, 
together with a partial slip BC, characterized by a non 
vanishing slip length (here b ~ llcr). An important point 
here is that the first layer of micro-ions now contributes a 
large amount to the global streaming current in contrast 
to the wetting emphasized in the inset of Fig. El 

The remobilization of the Stern layer therefore adds on 
to the slippage effect to increase the C, potential. 

We now summarize our results in Fig. |3and plot the 
C potential [deduced from the measure of the charge cur- 
rent supplemented with Eq. |j2Jl] as a function of the 
Debye screening factor in the wetting and non-wetting 
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FIG. 2: Measured Poiseuille velocity profile (solid line) in the 
wetting case [cps = !)• Dashed line : hydrodynamic pre- 
diction using a no-slip BC at the 'plane of shear' located at 
Zs (indicated by the arrow). To emphasize the existence of 
an immobile Stern layer, we also indicate the charge density 
profile pc{z) = e{p+{z) — p-{z)) (dotted line), with arbitrary 
units. The position of the wall (defined as that of the centers 
of the last layer of wall atoms) is at Zm = —10.8a. Inset: Re- 
sults for the non- wetting case {cfs = 0.5). Solid line : velocity 
profile measured in the simulation (shown on the same scale 
as in the main graph); Dashed line : hydrodynamic prediction 
with a partial slip BC, with a slip length b — 11a; Dashed- 
dot line : hydrodynamic prediction with a no-slip BC; Dotted 
line : charge density profile (arbitrary units). 
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cases. In this figure the C potential is normahzed by the 
bare surfacepotential Vq , obtained from the analytic PB 
expression [jj] (see Fig. Moreover, the slip length in 
the non- wetting case has been measured to barely depend 
on the screening factor, indicating that the micro-ions do 
not affect the fluid-solid friction. The overall conclusion 
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FIG. 3: The symbols show the potential measured in MD as 
a function the screening factor k£b, for the wetting (bottom) 
and non- wetting substrate (top). The ( potential is normal- 
ized by the bare surface potential Vb obtained from the PB 
expression at a given k and surface charge (see the discussion 
on the inset of Fig. 0. For the wetting case (bottom), the 
dashed line is the PB electrostatic potential V{zs) where the 
'plane of shear' position Zs does not vary significantly with 
salt. For the non- wetting case (top), the dash-dotted line 
corresponds to the slip prediction, Eq. Q (with b — 11a). 

from Fig. I^lis that non-wettability strongly amplifies the 
electrokinetic effects : the ratio between the ( potential 
and the surface potential is much larger in the hydropho- 
bic case as compared to the hydrophilic case. More pre- 
cisely, in the wetting case the ( potential is fixed by the 
electric properties of the surface, and coincides with the 
electric potential at the 'plane of shear', ~ V{zs), as 
is usually assumed 0. This is demonstrated in Fig. 13 
where the simulation points for ^ are compared to the PB 
estimate for the electric potential, Vpb{zs), showing an 
overall very good agreement. Conversely, the ( potential 
in the non-wetting case is dominated by the slip effect and 
the immobile Stern layer is completely absent. The ef- 
fect of slip can be accounted for by considering the partial 
slip BC in the electrokinetic current le = J dS pc(,z)v{z) 
with pc{z) the charge density and v{z) the velocity pro- 
file characterized by a slip length h. Within linearized 
PB description, and for the present planar geometry, the 
result for the current 7,., Eq. (|3Jl, may then be written 
le = ^^(1 + iib)fo In the non- wetting case, this 
leads to 

C^Vo{l + Kb) (4) 

with Vq the bare potential of the surface This ex- 

pression is successfully compared in Fig. ^to simulation 
results. 



Finally we quote that we have tested a different elec- 
trokinetic geometry, corresponding to the more common 
electro-osmotic situation : flow response to the appli- 
cation of an electric field. The ^ potentials measured in 
this geometry (not shown here) are in full qualitative and 
quantitative agreement with the results in the present 
(streaming current) geometry, both on the wetting and 
non-wetting surfaces. 

To conclude, we have shown using MD simulations that 
the notion of C potential, the cornerstone in the descrip- 
tion of EDL dynamics, encompasses different physical 
mechanisms, depending on the wettability of the charged 
substrate. In the wetting situation, the C potential can 
be indeed directly related to surface charge properties, 
confirming hereby the traditional Stern layer picture . 
In contrast, for non- wetting substrates, electric and slip 
effects are strongly intricated, leading for the ( potential 
(Eq. 0)) to an amplification ratio of 1 + nb. Practically, 
this ratio can take large values, even for moderate slip 
lengths on the order of a few nanometers, since the Debye 
length usually lies in the nanometer range (to be spe- 
cific, slip lengths of the order of tens of nanometers have 
been reported for hydrophobic silanized glass surfaces 
[i(tI | ) . It would be moreover highly desirable to generalize 
Eq. (01) obtained for the planar case to other geometries 
(eg. spherical) where a new length (radius R) will com- 
pete with the debye (k~^) and slip (6) lengths, Eq. Q 
yielding the limiting behavior for b/R, k^^/R — > 0. This 
work also points to the difference between dynamical and 
static properties of the EDL, respectively characterized 
by the Zcta and bare surface potential, C ^^'i ^o- In- 
dependent static and dynamic measurements, using e.g. 
AFM or SFA, should therefore allow to probe the differ- 
ent mechanisms underlined in this study. Another route, 
with potential applications in micro-fiuidics, consists in 
polarizing the solid substrate, i.e. "imposing" Vq, which 
should lead to large electrokinetic effects on non- wetting 
substrates. 
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